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Abstract 

This paper investigates general properties of codes with the rank metric. First, we investigate asymp- 
totic packing properties of rank metric codes. Then, we study sphere covering properties of rank metric 
codes, derive bounds on their parameters, and investigate their asymptotic covering properties. Finally, 
we establish several identities that relate the rank weight distribution of a linear code to that of its dual 
code. One of our identities is the counterpart of the MacWilliams identity for the Hamming metric, and 
it has a different form from the identity by Delsarte. 

I. Introduction 

Although the rank has long been known to be a metric implicitly and explicitly (see, for example, 
[1]), the rank metric was tirst considered for error control codes (ECCs) by Delsarte [2]. The potential 
applications of rank metric codes to wireless communications [3], [4], public-key cryptosystems [5], and 
storage equipments [6], [7] have motivated a steady stream of works [6]-[18], described below, that focus 
on their properties. 

The majority of previous works focus on rank distance properties, code construction, and efficient 
decoding of rank metric codes, and the seminal works in [2], [7], [8] have made significant contribution to 
these topics. Independently in [2], [7], [8], a Singleton bound (up to some variations) on the minimum rank 
distance of codes was established, and a class of codes achieving the bound with equality was constructed. 
We refer to this class of codes as Gabidulin codes henceforth. In [2], [8], analytical expressions to compute 
the weight distribution of linear codes achieving the Singleton bound with equality were also derived. 
In [6], it was shown that Gabidulin codes are optimal for correcting crisscross eiTors (referred to as 
lattice-pattern errors in [6]). In [7], it was shown that Gabidulin codes are also optimal in the sense 
of a Singleton bound in crisscross weight, a metric considered in [7], [11], [19] for crisscross errors. 
Decoding algorithms that parallel the extended Euclidean algorithm and the Peterson-Gomstein-Zierler 
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algorithm were introduced for Gabidulin codes in [8] and [7], respectively. In [2], the counterpart of 
the Mac Williams identity, which relates the rank distance enumerator of a code to that of its dual code, 
was established using association schemes. Following the works in [2], [7], [8], the construction in [8] 
was extended in [13] and the properties of subspace subcodes and subfield subcodes were considered 
in [14], [20]; the counterparts of the Welch-Berlekamp algorithm and the Berlekamp-Massey algorithm 
were considered in [21] and [22] respectively for Gabidulin codes; the error performance of Gabidulin 
codes was investigated in [11], [16], [17]. 

Some previous works focus on the packing and covering properties of rank metric codes. Both packing 
and covering properties are significant for ECCs, and packing and covering radii are basic geometric 
parameters of a code, important in several respects [23]. For instance, the covering radius can be viewed 
as a measure of performance: if the code is used for error correction, then the covering radius is the 
maximum weight of a correctable error vector [24]; if the code is used for data compression, then the 
covering radius is a measure of the maximum distortion [24]. The Hamming packing and covering radii of 
ECCs have been extensively studied (see, for example, [25]-[27]), whereas the rank packing and covering 
radii have received relatively little attention. It was shown that nontrivial perfect rank metric codes do 
not exist in [9], [10], [15]. In [12], a sphere covering bound for rank metric codes was introduced. 
Generalizing the concept of rank covering radius, the multi-covering radii of codes with the rank metric 
were defined in [28]. Bounds on the volume of balls with rank radii were also derived [18]. 

In this paper, we investigate packing, covering, and rank distance properties of rank metric codes. The 
main contributions of this paper are: 

• In Section |llll we establish further properties of elementary linear subspaces (ELS's) [17], and 
investigate properties of balls with rank radii. In particular, we derive both upper and lower bounds 
on the volume of balls with given rank radii that are tighter than their respective counterpart in [18]. 
These technical results are used later in our investigation of the sphere covering properties of rank 
metric codes. 

• In Section UVl we study the packing properties of rank metric codes, and also derive the asymptotic 
maximum code rate for a code with given relative minimum rank distance. 

• In Section lYl we study the covering properties of rank metric codes, and derive both upper and lower 
bounds on the minimal cardinality of a code with given length and rank covering radius. Our new 
bounds are tighter than the bounds introduced in [12]. Using the sphere covering bound, we also 
establish additional sphere covering properties for linear rank metric codes, and prove that some 
classes of rank metric codes have maximal covering radius. Finally, we establish the asymptotic 
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minimum code rate for a code with given relative covering radius. 

• In Section |VlJ we study the rank weight properties of linear codes. We show that, similar to the 
MacWilliams identity for the Hamming metric, the rank weight distribution of any linear code 
can be expressed as an analytical expression of that of its dual code. It is also remarkable that our 
MacWilliams identity for the rank metric has a similar form to that for the Hamming metric. Despite 
the similarity, our new identity is proved using a different approach based on linear spaces. The 
intermediate results obtained using our approach offer interesting insight. We also derive identities 
that relate moments of the rank weight distribution of a Unear code to those of its dual code. 

We provide the following remarks on our results: 

1) The concept of elementary linear subspace was introduced in our previous work [17]. It has similar 
properties to those of a set of coordinates, and as such has served as a useful tool in our derivation 
of properties of the rank metric (see Section HID), covering properties of rank metric codes in general 
(see Section IVj), and Gabidulin codes (see [17]). Although our results may be derived without the 
concept of ELS, we have adopted it in this paper since it enables readers to easily relate our 
approach and results to their counterparts for Hamming metric codes. 

2) Both the matrix form [2], [7] and the vector form [8] for rank metric codes have been considered 
in the literature. Following [8], in this paper the vector form over GF{q^) is used for rank metric 
codes although their rank weight is defined by their corresponding mxn code matrices over GF(g) 
[8]. The vector form is chosen in this paper since our results and their derivations for rank metric 
codes can be related to their counterparts for Hamming metric codes. 

3) In [2], the MacWilliams identity is given between the rank distance enumerator sequences of two 
dual codes using the generalized Krawtchouk polynomials. Based on a different proof, we establish 
the same identity for linear rank metric codes, although our identity is expressed using different 
parameters which are shown to be the generalized Krawtchouk polynomials as well. We also present 
this identity in weight enumerator polynomial form (cf. Theorem |2]). In their polynomial forms, 
the MacWilliams identities for both the rank and the Hamming metrics are similar to each other. 
Furthermore, the polynomial form allows us to derive further identities (cf. Propositions [18] and 
[T9b between the rank weight distribution of linear dual codes. 

The rest of the paper is organized as follows. Section HIl gives a brief review of necessary background 
to keep this paper self-contained. In Section [Till we derive some further properties of ELS's and balls 
of rank radii. In Sections |IV] and |Vl we investigate the packing and covering properties respectively of 
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rank metric codes. Section |Vl] presents our results on the rank weight distribution of rank metric codes. 

II. Preliminaries 

A. Rank metric and elementary linear subspaces 

Consider an n-dimensional vector x = (xo, xi, . . . , G GF((7'")". The field GF{q^) may be 

viewed as an m-dimensional vector space over GF(g). The rank weight of x, denoted as rk(x), is 
defined to be the maximum number of coordinates in x that are linearly independent over GF(q) [8]. 
Note that all ranks are with respect to GF(g) unless otherwise specified in this paper. The coordinates 
of X thus span a linear subspace of GF(g'"), denoted as S(x), with dimension equal to rk(x). For all 

det 

X, y G GF(g™)", it is easily verified that dR(x, y) = rk(x — y) is a metric over GF(g™)", referred to 
as the rank metric henceforth [8]. The minimum rank distance of a code C, denoted as df^{C), is simply 
the minimum rank distance over all possible pairs of distinct codewords. When there is no ambiguity 
about C, we denote the minimum rank distance as d^. 

In [17], we introduced the concept of elementary linear subspace (ELS). If there exists a basis set B 
of vectors in GF(g)" for a linear subspace V C G¥{q"^)^, we say V is an elementary linear subspace 
and B is an elementary basis of V. We denote the set of all ELS's of GF(g™)" with dimension v as 
Ev{q"^,n). An ELS has properties similar to those for a set of coordinates [17], and they are summarized 
as follows. A vector has rank < r if and only if it belongs to some ELS with dimension r. For any 
V G Ey{q'^,n), there exists V G £'„_t,(g™', n) such that V® V = GF(g'")", where © denotes the direct 
sum of two subspaces. For any vector x G GF(g™)", we denote the projection of x on V along V as 
xy, and we remark that x = xy + x-p. 

B. The Singleton bounds 

It can be shown that dg_ < du [8], where d^ is the minimum Hamming distance of the same code. 
Due to the Singleton bound for block codes, the minimum rank distance of an (n, k) block code over 
GF{q"^) thus satisfies 

4<n-A: + l. (1) 
An alternative bound on the minimum rank distance is also given in [29]: 

TTl 

dR< —(n-k) + l. (2) 
n 

For n < m, the bound in ^ is tighter than that in Q. When n > m the bound in Q is tighter. 
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When n < m, a class of codes satisfying ([T]) with equality was first proposed in [8] and then generalized 
in [13]. Let g = {go,gi, ■ ■ ■ ,9n-i) be linearly independent elements of GF(g™), then the code defined 
by the generator matrix 
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9n-l I 



(3) 



where [i] = q°-^ with a being an integer prime to m, is called a generalized Gabidulin code generated by 
g = ((Jo, 5i, . . . , 9n-i)\ it has dimension k and minimum rank distance = n — A; + 1 [13]. 

A class of codes satisfying Q with equality was proposed in [29]. It consists of cartesian products of a 
generalized Gabidulin code with length n = m. Let ^ be an (m, A;, c^r = m — + generalized Gabidulin 
code over GF(g'"), and let ^' =^ C/ x . . . x be the code obtained by / cartesian products of Q. Thus G'' 
is a code over GF((7™) with length ml, dimension kl, and minimum rank distance dj^ = m — k + 1 [29]. 



C. Covering radius and excess 

The covering radius p of a code C with length n over GF(g™) is defined to be the smallest integer p 
such that all vectors in the space GF{q"^)"- are within distance p of some codeword of C [27]. It is the 
maximal distance from any vector in GF((7™)" to the code C. That is, p = maXxgGF(g"')"{c^(x, C)}. 
Also, if C C C", then the covering radius of C is no more than the minimum distance of C". Finally, 
a code C with length n and minimum distance d is called a maximal code if there does not exist any 
code C with same length and minimum rank distance such that C C C". A maximal code has covering 
radius p < d — 1. 

Van Wee [30], [31] derived several bounds on codes with Hamming covering radii based on the excess 
of a code, which is determined by the number of codewords covering the same vectors. Below are some 
key definitions and results in [30], [31]. For all V C GF((/'")" and a code C with covering radius p, the 
excess on y by C is defined to be 

i?c(l^) = (4) 

CGC 

where Bp{c) denotes a ball centered at c with Hamming radius p. The excess on GF{q"^)"- by C is given 
by £^(7(GF(g™)") = [C] -l/f (q'™, n) -q'™", where Vpiq"", n) denotes the volume of a ball with Hamming 
radius p. Also, if {Wi] is a family of disjoint subsets of GF(g™)", then £^c(Ui^i) = T^i^ciWi). 

def 

Suppose Z = {z e GF(g™)"|£'c({z}) > 0} [30], i.e., Z is the set of vectors covered by at least two 
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codewords in C. Note that z G Z if and only if |S^(z) n C| > 2. It can be shown that \Z\ < Ec{Z) = 
Sc(GF(g™)") = \C\ ■ V^{q'^,n) - q""". 

Although the above definitions and properties were developed for the Hamming metric, they are in 
fact independent of the underlying metric and thus are applicable to the rank metric as well. 

D. Notations 

In order to simplify notations, we shall occasionally denote the vector space GF(g'")" as F. We denote 
the number of vectors of rank u (0 < u < min{m,n}) in GF((7"^)" as Nu{q"^,n). It can be shown that 
Nu{q'^,n) = [I] a(m, u) [8], where a(m, 0) = 1 and a(m, u) = lYi=o il"" - Q') for ^ > 1. The term 
is often referred to as a Gaussian polynomial [32], defined as ["] =^ a(n, u)/a{u, u). Note that ["] is the 
number of n-dimensional linear subspaces of GF(g)". We refer to all vectors in GF(g™)" within rank 
distance r of x G GF{q"^)"' as a ball of rank radius r centered at x, and denote it as i3r(x). Its volume, 
which does not depend on x, is denoted as Vr{q''^,n) = Ylu=o ^^il"^ jn). We also define /?(m, 0) =' 1 
and P{m, u) = nr=o ["Tl for ^ > 1' which are used in Section |Vll These terms are closely related to 
Gaussian polynomials: /3(m, u) = [™]/?(m, n) and /3(m + u, m + n) = n). 

III. Technical results 

A. Further properties of ELS' s 

Lemma 1: Any vector x € GF(g™)" with rank r belongs to a unique ELS V € Er{q"^,n). 

Proof: The existence of V € Er{q'^, n) has been proved in [17]. Thus we only prove the uniqueness 
of V, with elementary basis {vj}[Z^(}. Suppose x also belongs to W, where W € £'r(g'",n) has an 
elementary basis {^jYjZ^- Therefore, x = X]i=o — Yl^j=o^j'^j^ where ai,bj € GF(g™) for < 
hj < r — 1. By definition, we have ©(x) = &{ao, . . . , a,.-i) = &{bo, . . . , br-i), therefore 6j's can be 
expressed as linear combinations of a^'s, i.e., bj = Y^^Iq Cj^ai where Cj^i € GF(g). Hence 

r— 1 r— 1 r— 1 r— 1 

X = ^ bjV^j = CjiOiVirj = ^ OjUj, (5) 

j=0 j=0 i=0 j=0 

where Uj = X]j=o '^jA'^j ^ GF(g)". Now consider X, the matrix obtained by expanding the coordinates 
of X with respect to the basis {aj}™Q^. For < i < r — 1, the i-th row of X is given by the vector Vj 
by definition and by the vector Uj from Eq. Therefore = Uj € W, and hence V C W. However, 
dim(V) = dim(W), and thus V = W. ■ 
Lemma [T] shows that an ELS is analogous to a subset of coordinates since a vector x with Hamming 
weight r belongs to a unique subset of r coordinates, often referred to as the support of x. 
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In [17], it was shown that an ELS always has a complementary elementary linear subspace. The 
following lemma enumerates such complementary ELS's. 

Lemma 2: Suppose V € Ey{q'^,n) and ^ C V is an ELS with dimension a, then there are q"-^'"~'^^ 
ELS's B such that A(B B = V. Furthermore, there are [^] such ordered pairs {A, B). 

Proof: First, remark that dim(;B) = v — a. The total number of sets of f — a linearly independent 
vectors over GF{q) in is given by = {q^ — q"-){q'" — q'^'^'^) ■ ■ ■ {q^ — q^~^) = q'^^'"~°'^a{v — a, v — a). 
Note that each set of linearly independent vectors over G¥{q) constitutes an elementary basis set. Thus, 
the number of possible B is given by N divided by a{v — a,v — a), the number of elementary basis sets for 
each B. Therefore, once A is fixed, there are q'^i'"-"-) choices for B. Since the number of a-dimensional 
subspaces ^ in V is [^] , the total number of ordered pairs {A, B) is hence q^-""") [^] . ■ 

Puncturing a vector with full Hamming weight results in another vector with full Hamming weight. 
Lemma [3] below shows that the situation for vectors with full rank is similar. 

Lemma 3: Suppose V € E^{q^"', n) and u € V has rank v, then rk(u_4) = a and rk(ug) = v — a fox 
any A € n) and B G E^^aiq"", n) such that A(S)B = V. 

Proof: First, u_4 E A and hence rk(u^) < a by [17, Proposition 2]; similarly, rk(ug) < v — a. 
Now suppose rk(u_4) < a or rk(ug) < v — a, then v = rk(u) < rk(u^) + rk(u0) <a + v — a = v. ■ 

It was shown in [17] that the projection u_4 of a vector u on an ELS A depends on both A and its 
complement B. The following lemma further clarifies the relationship: changing B always modifies u^, 
provided that u has full rank. 

Lemma 4: Suppose V € E^{q™-,n) and u G V has rank v. For any A G Ea{q"^,n) and B G 
Ey^a{Q"^,n) such that A (B B = V, define the functions f^i{A,B) = u_4 and g-^{A,B) = ug. Then 
both fu and axe injective. 

Proof: Consider another pair [A', B') with dimensions a and v — a respectively. Suppose A' ^ A, 
then u_4' 7^ u_4. Otherwise u_4 belongs to two distinct ELS's with dimension a, which contradicts 
Lemma [T] Hence u_4' ^ u_4 and ug' = u — u_4' 7^ u — u_4 = ug. The argument is similar if B' ^ B. ■ 

B. Properties of balls with rank radii 
Lemma 5: For < r < min{n, m}, 

where a{q) J2T=i k( ^-i) ^ decreasing function of q satisfying a{q) < 2 for g > 2 [17]. 
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Proof: The upper bound in Q was derived in [17, Lemma 13], and it suffices to prove the lower 
bound. Without loss of generality, we assume that the center of the ball is 0. For any x e G¥{q^Y , 
we associate one subspace T of G¥{q'^) such that dim(T) = r and S(x) C T. We consider the vectors 
y G GF(g'")"~'^ such that S(y) C T. There are q""' choices for x and, for a given x, q^^^~^^ choices 
for y. Thus the total number of vectors z = (x, y) e GF((7™)" is Since S(z) C T, we have 

rk(z) < r. Thus, K(g™,n) > qr(m+n-T) _ , 

We remark that both bounds in Q are tighter than their respective counterparts in [18, Proposition 1]. 
More importantly, the two bounds in Q differ only by a factor of g'^^'?), and thus they not only provide 
a good approximation of Vr{q^,'n), but also accurately describe the asymptotic behavior of Vr{q'^,n). 

The diameter of a set is defined to be the maximum distance between any pair of elements in the set 
[25, p. 172]. For a binary vector space GF(2)" and a given diameter 2r < n, Kleitman [33] proved that 
balls with Hamming radius r maximize the cardinality of a set with a given diameter. However, when 
the underlying field for the vector space is not GF(2), the result is not necessarily valid [27, p. 40]. We 
show below that balls with rank radii do not necessarily maximize the cardinality of a set with a given 
diameter. 

Proposition 1: For 3 < n < m and 2 < 2r < n, there exists S C GF((7™)" with diameter 2r such 

that \S\ > 

Proof: The set S = {(xq, . . . , G GF(g'")"|x2r = ••• = = 0} has diameter 2r 

and cardinality g^™'". For r = 1, we have Vi{q^,n) = 1 + ^'^ ~{q-i) ^ 1^^- r > 2, we have 
Vriq"", n) < qr{n+m)-(r^-a{q)) LemmaS Since > 2 > a{q), we obtain Vr{q'^,n) < q^i^+^) < \S\. 

m 

The intersection of balls with Hamming radii has been studied in [27, Chapter 2], and below we 
investigate the intersection of balls with rank radii. 

Lemma 6: If < r, s < n and ci, C2 G GF((7™)", then \Br{ci) n Bs{c2)\ depends on ci and C2 only 
through dR(ci,C2). 

Proof: First, without loss of generality, we assume ci = 0, and we denote rk(c2) = e. We can 
express C2 as C2 = uB, where u = (uq, . . . , Ug-i, 0, . . . , 0) G GF(g'")" has rank e and B G GF((7)"^"' 
has full rank. For any x G i?r(0)ni?s(u) we have rk(xB) = rk(x) < r and rk(xB— C2) = rk(x— u) < s. 
Thus there is a bijection between 5^(0) n -Bs(uB) and Br{0) n -Bs(u). Hence \Br(0) n 5s(uB)| = 
l-Br(O) n -^^(u)!, that is, [-6^(0) n ^^(uB)! does not depend on B. 

Since |i?r(0) PI i?s(uB)| is independent of B, we assume B = Inxn without loss of generality 
henceforth. The nonzero coordinates of u all belong to a basis set {ui}^^ of GF(g'"). Let x = 
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(xo, • • • ,Xn-i) € Br{0) n Bs{u), then we can express Xj as xj = cnjUi with ajj € GF(g) for 

< i < n — 1. Suppose v = {vq, . . . , fe-i, 0, . . . , 0) G GF(g'")" also has rank e, then the nonzero 
coordinates of v all belong to a basis set {wj}[!lo^ GF(g'"). We define x = {xq, . . . , G GF(g™)'^ 

such that Xj = J2Y=o^ CLijVi for < j < n — 1. We remark that rk(x) = rk(x) < r and rk(x — v) = 
rk(x— u) < s. Thus there is a bijection between Br{0)nBs{yr) and Br{0)nBs{u). Hence \Br{0)nBs{u)\ 
depends on the vector u only through its rank e. ■ 
Proposition 2: If < r,s < n, ci, C2, c'^^, C2 G GF(g™)" and ^5(01,02) > df,{c[,C2), then 

|B,(ci)n5,(c2)I < |5,(c;)nB,(c'2)|. (7) 

Proof: It suffices to prove (|7]) when ^^(01,02) = df^{c[,C2) + 1 = e + 1. By Lemma [6j we 
can assume without loss of generality that ci = c[ = 0, C2 = (0, ci, . . . , Ce, 0, . . . , 0) and C2 = 
(co, ci, . . . , Ce, 0, . . . , 0), where cq, . . . , Ce G GF(q'™) are linearly independent. 

We will show that an injective mapping (p from Br{ci)f]Bs{c2) to Br{c[)nBsic2) can be constructed. 
We consider vectors z = (zq, zi, . . . , Zn-i) G -^^(ci) n Bs{c2). We thus have rk(z) < r and rk(u) < s, 
where u = {uq, ui, . . . , tin-i) = z — C2 = {zq — cq^zi — ci, . . . , Zn-i)- We also define z = (zi, . . . , 2;„_i) 
and u = (ui, . . . , Un-i)- We consider three cases for the mapping (j), depending on z and u. 

• Case I: rk(u) < s — 1. In this case, (/)(z) z. We remark that rk(z — C2) < rk(u) + 1 < s and 
hence (/)(z) G -Br(ci) n -Bs(c2). 

_ _ def 

• Case II: rk(u) = s and rk(z) < r — 1. In this case, </)(z) = {zq — cq, zi, . . . , Zn^i). We have 
rk(0(z)) < rk(z) + l < r and rk ((/)(z) - 03) = rk(z - C2) < s, and hence (^(z) G Br{c'^)<r\Bs{c'2). 

• Case III: rk(u) = s and rk(z) = r. Since rk(u) = s, we have zq — cq ^ Similarly, since 
rk(z) = r, we have zq G S(z). Denote dim(©(u, z)) as d (d > s). For d > s, let aoi • • • > Orf-i be a 
basis of ©(u, z) such that oq, . . . , as-i £ and q^, . . . , Ud-i G S(z). Note that cq G G(u, z), 
and may therefore be uniquely expressed as cq = c„ + c^, where c„ G S(ao, • • • ^ 

_ _ def 

and Cz G ©(a^, . . . , a^-i) ^ ©(z)- If d = s, then = G ©(z). In this case, 0(z) = {zq — 
Cz,zi, . . . , Zn-i). Remark that zq — Cz G ©(z) and hence rk((?!<(z)) = r. Also, zq — Cz = zq — cq + Cu G 
©(u) and hence rk(</)(z) — C2) = s. Therefore </)(z) G Br{c[) n i?s(c2). 
We now verify that the mapping </> is injective. Suppose there exists z' such that (/){z') = <j){z). Since 
0(z) only modifies the first coordinate of z, the last n — 1 coordinates of z and z' are equal and so are 
the last n — 1 coordinates of z — C2 and z' — C2. Hence z and z' belong to the same case. It can be easily 
verified that for each case above, 4> is injective. Hence 0(z') = (/)(z) implies that z' = z. Therefore (/> is 
injective, and \Br{ci) n -65(02) [ < l-Br(c'i) n -65(02) |. ■ 
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Corollary 1: If <r,s <n, ci,C2,c'^,C2 € GF(g™)" and (iR(ci,C2) > ^^(c']^, Cg), then 

\Brici) U B,(C2)I > \Bric[) U 5.(4)1. (8) 
Proof: The result follows from |Sr(ci) U5,(c2)| = ^^(9™, + ^.(g"", - |Sr(ci) n Ss(c2)|. ■ 
We now quantify the volume of the intersection of two balls with rank radii for some special cases, 
which will be used in Section IV-BI 

Proposition 3: If ci, C2 G GF(g™)" and 4(ci, C2) = r, then |-Br(ci) n Si(c2)| = 1 + (g™ - q"^) [[] + 

Proof: The claim holds for r = m trivially, and we assume r < m henceforth. By Lemma |6l 
assume C2 = and hence rk(ci) = r without loss of generality. By Lemma [T] the vector ci belongs 
to a unique ELS V G Er{q^,n). First of all, it is easy to check that y = G -Br(ci) n -Bi(O). We 
consider a nonzero vector y G i?i(0) with rank 1. Firstly, if y G V, then ci — y G V. We hence have 
rk(ci — y) < r and y G Bj.{ci). Note that there are {q^ — 1) such vectors. Secondly, if y ^ V and 
®(y) C <S(ci), then 6(ci — y) C ©(ci). We hence have rk(ci — y) < r and y G Br{ci). Note that 
there are (g*" — !)(["] — [^] ) such vectors. Finally, suppose y ^ V and S(y) ^ S(ci). Denote the linearly 
independent coordinates of ci as ao, • • • ,ar-i and a nonzero coordinate of y as ^ S(ci), where 
{oj}™ is a basis set of GF((7™). Then the matrix Ci — Y obtained by expanding the coordinates of 
ci — y according to the basis {oj} has row rank r + 1. Therefore rk(ci — y) = r + 1, and y ^ Br{ci). 

■ 

Proposition 4: If ci,C2 G GF(g'")" and 4(ci,C2) = r, then \Bs{ci) f^ Br-s{c2)\ = g''^''"''^ Q for 
< s < r. 

Proof: By Lemma [6l we can assume that ci = 0, and hence rk(c2) = r. By Lemma [T] C2 belongs 
to a unique ELS V G Er{q"^,n). We first prove that all vectors y G Bs{0) n i3r-s(c2) are in V. Let 
y = yv + yw, where W G En-riq"", n) such that V W = GF(g'")". We have yy + (c2 - y)v = C2, 
with rk(yv) < rk(y) < s and rk((c2 — y)v) < i'k(c2 — y) < r — s. Therefore, rk(yv) = rk(y) = s, 
rk((c2 - y)v) = rk(c2 - y) = r - s, and 6(yv) n 6((c2 - y)v) = {0}. Since rk(yv) = rk(y), we 
have 6(y>v) ^ ®(yv); and similarly ©((c2 — y)>v) ^ ©((^2 — y)v)- Altogether, we obtain 6(y>v) H 
S((c2 - y)>v) = {0}. However, yw + (c2 - y)w = 0, and hence y^ = (c2 - y)w = 0. Therefore, 
y G V. 

We now prove that y is necessarily the projection of C2 onto some ELS ^ of V. If y G V satisfies 
rk(y) = s and rk(c2 — y) = r — s, then y belongs to some ELS A and C2 — y G S such that A® B = V. 
We hence have y = C2,a and C2 — y = C2,e- 
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On the other hand, for any A G Es{q"^,n) and B G Er-siQ"^,^) such that A (B B = V, C2,A is a 
vector of rank s with distance r — s from C2 by Lemma [3l By Lemma IH all the C2.a vectors are distinct. 



Fig. 1. Illustration of Proposition [4] 

As shown in Figure [T] only the outmost layers of two balls of radii s and r — s intersect when the 
distance between the two centers is r. Proposition |4] quantifies the volume of the intersection in Figure [1] 

The problem of the intersection of three balls with rank radii is more complicated since the volume 
of the intersection of three balls with rank radii is not completely determined by the pairwise distances 
between the centers. We give a simple example to illustrate this point: consider GF(2^)'^ and the vectors 
ci = c'^ = (0, 0, 0), C2 = C2 = (1, a, 0), C3 = (a, 0, 1), and Cg = (a, a + 1, 0), where a is a primitive 
element of the field. It can be verified that d^{c,\^C2) = (iR(c2,C3) = dn{c^,ci) = 2 and dR{c[,C2) = 
d^(c'2,C3) = 4(0^5, c'l) = 2. However, Si(ci) n ^1(02) n -61(03) = {{a + 1,0,0)}, whereas Bi{c[) n 
^1(02) n ^1(03) = {(1,0,0), (0,a + 1,0), (a,a,0)}. We remark that this is similar to the problem of 
the intersection of three balls with Hamming radii discussed in [27, p. 58], provided that the underlying 
field GF{q'^) is not GF(2). 

IV. Packing properties of rank metric codes 

Combining ([U and (|2]l and generalizing slightly to account for nonlinear codes, we can show that the 
cardinality K of a code C over GF(g'") with length n and minimum rank distance dg_ satisfies 



In this paper, we call the bound in ^ the Singleton bounqj for codes with the rank metric, and refer to 
codes that attain the Singleton bound as maximum rank distance (MRD) codes. 

'The Singleton bound in [7] has a different form since array codes are defined over base fields. 



There are thus as many vectors y as ordered pairs {A,B). By Lemma |2l there are g'^^'" [^] such pairs, 
and hence vectors y. ■ 





(9) 
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For any given parameter set n, m, and ^r, explicit construction for MRD codes exists. For n < m and 
dR < n, generalized GabiduUn codes can be constructed. For n > m and < m, an MRD code can 
be constructed by transposing a generalized GabiduUn code of length m and minimum rank distance c?r 
over GF(g"), and this code is not necessarily linear over GF(g"^). Although maximum distance separable 
(MDS) codes, which attain the Singleton bound for the Hamming metric, exist only for limited block 
length over any given field, MRD codes can be constructed for any block length n and minimum rank 
distance over arbitrary fields GF(g"*). This has significant impact on the packing properties of rank 
metric codes as explained below. 

The sphere packing problem we consider is as follows: given a finite field GF(g™), length n, and 
radius r, what is the maximum number of non-intersecting balls with radius r that can be packed into 
GF(g"*)'^? The sphere packing problem is equivalent to finding the maximum cardinality A{q'^, n, d) of 
a code over GF(g"*) with length n and minimum distance c? > 2r + 1: the spheres of radius r centered 
at the codewords of such a code do not intersect one another. Furthermore, when these non-intersecting 
spheres centered at all codewords cover the whole space, the code is called a perfect code. 

For the Hamming metric, although nontrivial perfect codes do exist, the optimal solution to the sphere 
packing problem is not known for all the parameter sets [25]. In contrast, for rank metric codes, although 
nontrivial perfect rank metric codes do not exist [9], [10], we show that MRD codes provide an optimal 
solution to the sphere packing problem for any set of parameters. For given n, m, and r, let us denote the 
maximum cardinality among rank metric codes over GF{q"^) with length n and minimum distance = 
2r + 1 as Af.{q"^, n, d^). For d^ > min{n, m}, ^R(g™, n, d^) = 1. For d^ < min{n, m}, ^R(g™', n, d^) = 
min |^'"("~'^'*+^), }. Note that the maximal cardinahty is achieved by MRD codes for all 
parameter sets. Hence, MRD codes admit the optimal solutions to the sphere packing problem for rank 
metric codes. 

The performance of Hamming metric codes of large block length can be studied in terms of asymptotic 
bounds on the relative minimum distance in the limit of infinite block length. In this section, we aim 
to derive the asymptotic form of AR{q^, n, d^) in the case where both block length and minimum rank 
distance go to infinity. However, this cannot be achieved for finite m since the minimum rank distance 
is no greater than m. Thus, we consider the case where lim„_»oo m ~ ^' where 6 is a constant. 

Define S == lim„^oo ^ and a{6) == lim„^oo sup '"Sg" .".W) ^ where a{6) represents the 
maximum possible code rate of a code which has relative minimum distance S as its length goes to 
infinity. We can thus determine the maximum possible code rate a{6) of a code. 

Proposition 5: For < 6 < min{l, b~^}, the existence of MRD codes for all parameter sets impUes 
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that 

a((5) = mill {1- (5,1- M}. (10) 

V. Covering properties of rank metric codes 

A. The sphere covering problem 

In this section, we are interested in the sphere covering problem for the rank metric. This problem can 
be stated as follows: given an extension field GF(g'"), length n, and radius p, we want to determine the 
minimum number of balls of rank radius p which cover GF(g'")" entirely. The sphere covering problem 
is equivalent to finding the minimum cardinality K^{q^ ,n, p) of a code over GF(g'") with length n and 
rank covering radius p. We remark that if C is a code over GF(g"^) with length n and covering radius p, 
then its transpose code is a code over GF(g") with length m and the same covering radius. Therefore, 
K^{q'^ ,n, p) = Kf^{q^ p), and without loss of generality we shall assume n < m henceforth in this 
section. 

We remark that Kg,{q"^,n,0) = g'"" and Kg,{q'^,n,n) = 1 for all m and n. Hence we assume 
< p < n throughout this section. Two bounds on Ky^{q^ ,n, p) can be easily derived. 

Proposition 6: For a code over GF(g'") with length n and covering radius < p < n, we have 

Trhr-^ +l<K,{q^,n,p)<q^(^-P\ (11) 
Proof: The lower bound is a straightforward generalization of the bound given in [12]. Note that 
the only codes with cardinality y are perfect codes. However, there are no nontrivial perfect codes 

for the rank metric [9]. Therefore, K^{q^,n^ p) > y . The upper bound follows from p < n — k for 
any (n, k) linear code [27], and hence any linear code with covering radius p has cardinality < q^^^~P) . 

■ 

We refer to the lower bound in ([TT] ) as the sphere covering bound. 

For a code over GY{q^) with length n and covering radius < p < n, we have K^{q^ p) < 
Kn{q"^, n, p), where Kii{q^,n, p) is the minimum cardinality of a (linear or nonlinear) code over GF((7™) 
with length n and Hamming covering radius p. This is because any code with Hamming covering radius 
p has rank covering radius < p. Since Kn{q'^ ,n, p) < q"^("'~p), this provides a tighter bound than the 
one given in Proposition [6] 

Lemma 7: For all m > and nonnegative n, n', p, and p', we have 

n + n', p + p') < K^iq"', n, p)K^{q"', n' , p'). (12) 
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In particular, we have 

K^{q^,n + l,p+l) < K^{q'^,n,p), (13) 

K,(g'-,n+l,p) < (14) 
Proof: First, (fT2l ) follows directly from [12, Proposition 4]. In particular, when {n',p') = (1, 1) and 
{n',p') = (1,0), we obtain ([T3] ) and ([141 ) respectively. ■ 



B. Lower bounds for the sphere covering problem 

We will derive two nontrivial lower bounds on K^{q"^ ,n, p). First, we adapt the bound given in [34, 
Theorem 1]. 

Proposition 7: For all g™, n, and < p < n, we have 

'q^^-A^{q^,n,2p+l)qP' ^ 



K,{q'^,n,p)> 



Vpiq"',n)-qp' 



2pl 



(15) 



provided that the denominator on the right hand side (RHS) is positive. 

Proof: Suppose C is a code over GF{q"^) with length n and rank covering radius p, and let Co be a 
maximal subcode of C with minimum rank distance d' > 2p+ 1. If d' > n, we choose Co to be a single 
codeword. Co thus covers \Co\Vp{q"^ , n) vectors. Define Ci = C\Co, (Ci is not empty, otherwise C 
would be a nontrivial perfect code) and for any Ci € Ci, let /(ci) denote the number of vectors covered 
by ci which are not covered by Co. Since Co is maximal, there exists at least one codeword co G Co 
such that (iR(co, ci) < 2p. We have /(ci) < Vp{q"^, n) — qP^ \^p'\ , where the equality corresponds to when 
there is only one such cq and d^{cQ,ci) = 2p by Proposition |2] In that case. Proposition |4] implies that 
there are qP^ [^^] vectors covered by both cq and ci. Thus, we have 

g"^" < n\Vp{q'^,n)+ /(ci) 

CiGCi 

2p 
P . 



< |Co|y,(g™,n) + (|C| - |Co|) Vpiq'^^n) - q 



\C\[Vp{q^,n)-qP 



■2p 
P . 



+ \Co\q'' 



2p 
P . 



We have \Cq\ < n, d') < ^(9"", n, 2p + 1), and the result follows. ■ 

We remark that A^{q'^,n,2p + 1) is if 2p + 1 < n, or 1 otherwise. Next, we obtain both 

sufficient and necessary conditions under which the bound is nontrivial, i.e., when the denominator on 

the RHS of ([15]) is positive. 

Lemma 8: The denominator on the RHS of ([TSl) is positive if p{m + n — 3p) > (T{q). Also, the 

denominator in (TSh is positive only if m + n > 3p. 
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Proof: We first prove the sufficient condition. We need to show Vp{q^, n) > qP^ [^^] . By Lemma |5l 
Vpiq'^^n) > gp("^+"-p). By [17, Lemma 1], we have [^p] < qp'+'^^i). Therefore, the denominator in ^ 
is positive if p{m + n — p) > 2/9^ + a{q). 

We now prove the necessary condition. Note that a{n,p) < q^P and a{2p,p) > g2p--T{g)^ where 
r(g) = logq (^^tt) [17, Lemma 2]. Now suppose p{m + n - 3p) < -T{q), then qP{r'^+n-3p)+r(q) ^ 
1. This implies ^^9^^™"^^ < 1, and hence ['^]q"'P < [^l^]q'''- By [17, Lemma 13], we obtain 
Vpiq"",?!) < [^]q'^P < qp' ^P]. Therefore, Vp{q'^,n) - qp" ^P] > only if p{m + n - 3p) > -T{q). 
Finally, < r(g) < 1 for q > 2 and hence p(m + n — 3p) > 0. ■ 

Before deriving the second nontrivial lower bound, we need the following adaptation of [31, Lemma 

def 

8]. Let C be a code with length n and rank covering radius p over GY{q'^). We define A = {x G 
GF((7'"r|4(x,C) =p}. 
Lemma 9: For x E A\Z and < p < n, we have 



(16) 



where 



def 



■p+il 



1 



Proof: Since x ^ Z, there is a unique cq G C such that (iR(x, cq) = p. By Proposition [3] we have 
\Bp{co)nBi{x)\ = l + Forany codeword ci € C satisfying 4(x, ci) = p+1, 

by Proposition |4] we have \Bp{ci)r]Bi (x) | = qP [''^^] . Finally, for all other codewords C2 G C at distance 

def 

> p + 1 from X, we have \Bp{c2) ni?i(x)| = 0. Denoting N = \{ci G C|(iR(x, ci) = p+ we obtain 



EciBiix)) 



^|Sp(c)n5i(x)| - 1^1 (x 

cGC 



qP) 
qP) 



+ NqP 



1 



if 



'p 




n 


^ mod ^(7'' 




1 




1 




1 



The proof is completed by realizing that {q^ — qP) ( [^] — ["])< 0, while Ec{Bi{-yi)) is a non-negative 
integer. ■ 



Proposition 8: If e > 0, then 

K^{q"',n,p)> 



Vp{q"^,n)-INp{q--,n) 



(17) 



where 5 = Vi {q"^, n) - qP-^ [?] - 1 + 2e. 
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The proof of Proposition [H provided in Appendix lAl uses the approach in the proof of [31, Theorem 6] 
and is based on the concept of excess reviewed in Section III-CI We remark that, unlike the bound given 
in Proposition |7J the bound in Proposition [8] is always applicable. The lower bounds in ( fT5] ) and dTT] ). 
when applicable, are at least as tight as the sphere covering bound given in (fTTb . 

C. Upper bounds for the sphere covering problem 

From the perspective of covering, the following lemma gives a characterization of MRD codes in terms 
of ELS's. 

Lemma 10: Let C be an (n, k) linear code over GF(g"^) (n < m). C is an MRD code if and only if 
C e V = GF(g™)" for all V G ^„-fc(g™, n). 

Proof: Suppose C is an {n, k,n — k + 1) MRD code. It is clear that C H V = {0} and hence 
C e V = GF(g'")" for all V G E^-kiq"", n). 

Conversely, suppose C© V = GF(g™)" for all V € En-k{(f^ ■, n). Then C does not contain any nonzero 
codeword of weight < n — k, and hence its minimum distance is n — k + I. ■ 

Let ao, ai, . . . , Om+p-i G GF(g™+'') be a basis set of GF(g"+^) over GF{q), and let /3o, • • • , Pm-i 
be a basis of GF{q"^) over GF(g). We define the linear mapping / between two vector spaces GF{q"^) 

del 

and Gm = S(ao, ai, . . . , am-i) by /(/?«) = Oj for < i < m — 1. This can be generalized to n- 
dimensional vectors, by applying / componentwise. We thus define / : GF{q"^)"- —>■ GF{q"^'^P)"' such 
that for any v = {vq, . . . , Vn-i), /(v) = {/{vq), . . . , f(vm-i))- This function / is a linear bijection from 
GF(q™)" to its image S;;. 
Lemma 11: For all r and any V G Er{q'^,n), /(V) C W, where W G Er{q"^~^^ ,n). 

Proof: We first show that / preserves the rank. Suppose u G GF(g™)". Let us denote the matrix 
formed after extending the coordinates of u with respect to the basis {Pi} as U. The extension of 

- (^\ 

/(u) with respect to the basis {ai} is given by U = . We thus have rk(U) = rk(U), and 

V / 

rk(/(u)) = rk(u). 

Let B = {hi} be a basis of V G Er{q"^,n) with vectors of rank one. Then for all i, /(bj) has rank 
one and {/(bj)} form a basis, and hence /(V) C W, where W G Er{q^~^P ,n) with {/(bj)} as a basis. 

■ 

Proposition 9: Let C be an {n,n — p, p + 1) MRD code with covering radius p. Then the code /(C) 
is a code of length n over GF((7'"'^'') with cardinality and covering radius p . 

Proof: The other parameters for the code are obvious, and it suffices to establish the covering 
radius. Let Tp be a subspace of GF{q"^^P) with dimension p such that 6m © = GF{q"^^P). Any 
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u € GF{q"^~^P)"' can be expressed as u = v + w, where v G &^ and w G T^. Hence rk(w) < p, and 
w E yV for some W G Ep{q"^^P, n). By Lemma [TOl we can express v as v = /(c + e) = /(c) + /(e), 
where c € C and e G V, such that / (V) C W. Eventually, we have u = /(c) + /(e) + w, where 
/(e) +w G W, and thus d(u, /(c)) < p. Thus /(C) has covering radius < p. Finally, it is easy to verify 
that the covering radius of / (C) is exactly p. ■ 
Corollary 2: We have 

KR(g'",n,p) < gmax{m-p,n}(n-p)_ (^g^ 

Proof: We can construct an (n,n — p) MRD code C over GF(g^) with covering radius p, where 
p = max{m — p,n}. By Proposition |9l /(C), where / maps GF(g^)" into a subset of GF(g'")", has 
covering radius p. Note that |/(C)| = |C| = ■ 
We can use the properties of Kf,{q"^,n, p) in Lemma |7] in order to obtain two tighter bounds when 
p > m — n. 

Proposition 10: Given fixed m, n, and p, for any n > I > and {ui, pi) for < i < / — 1 so that 
< rii < n, < Pi < Ui, and Ui + pi < m for all i, and X]!=o = and X]!=o Z'* ~ Z^' '^^^^ 

ifR(g™,n,«)< min /g"^("-p)-E.P.K-pO\ . (19) 

{(ni.pi): 0<i<Z-l} I J 

Proof: By Lemma |7J we have K^{q''"- ,n, p) < Y[iK!i{q"^,ni, pi) for all possible sequences {pi} 
and {nj}. For all i, we have Kf,{q^ ^m, pi) < qi(™^P»)("^P») by Corollary |2l and hence K^{q^ ,n, p) < 

It is clear that the upper bound in ([T9l ) is tighter than the upper bound in ([TTI ). It can also be shown 
that it is tighter than the bound in (fTSl ). 

The following upper bound is an adaptation of [27, Theorem 12.1.2]. 

Proposition 11: For any m, n < m, and p < n, there exists a code over GF(g'") of length n and 
covering radius p with cardinality 

1 



K^{q"',n,p) < 



+ 1. (20) 



l-logg,„„ (g'^"-yp(g'",n)) 
Our proof, given in Appendix iBl adopts the approach used to prove [27, Theorem 12.1.2]. 

Proposition 12: For all m, n < m, p < n, we have 

K^{q"',n,p) < [l + ln(y,(g",n))]. (21) 

Proof: Consider the square 0-1 matrix A of order g'"", where each row and each column represents 
a vector in GF(g™)". Set = 1 if and only if the sphere with rank radius p centered at vector i covers 
the vector j. There are thus exactly Vp{q'^,n) ones in each row and each column of A. Note that any 
qmn ^ submatrix C of A with no all-zeros rows represents a code with cardinality K and covering 
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radius p. Applying the Johnston-Stein-Lovasz theorem [27, Theorem 12.2.1] to A, we can find such a 
submatrix with K < y [g"*" + g™" IniVpiq"^, n))] . ■ 
The tightest bounds on K^{q^ p) known so far are given in Table U for q = 2, 2<m<l, 
2 < n < m, and 1 < p < 6. 

D. Covering properties of linear rank metric codes 

For a linear code with given covering radius, the sphere covering bound also imphes a lower bound 
on its dimension. 

Proposition 13: An {n,k) linear code over GF(g'") with rank covering radius p satisfies 

^_p_ P(--p)+^il) +i<k<n-p. (22) 
m J 

Proof: The upper bound directly follows the upper bound in ([TT]) . We now prove the lower bound. 
By the sphere covering bound, we have q"^'' > vJo^^Hij- However, by Lemma |5] we have Vp{q"^,n) < 

qPijn+n-p)+a(q) hcnCC > qmn—p{m+n-p]-a{q) _ _ 

We do not adapt the bounds in ([TSl l and (fTTl ) as their advantage over the lower bound in (l22l) is not 
significant. Table |ll] lists the values of the bounds in Proposition [13] on the dimension of a linear code 
with given covering radius for 4 < m < 8 and 4 < n < m. Note that only 1 < p < n — 1 are considered 
in Table JI] since, as shown below, the dimension of a linear code with given covering radius can be 
completely determined when p € {0, 1, n — 1, n}. 

Proposition 14: Let C be an {n,k) linear code over GF{q'^) (n < m) with rank covering radius p. 
Then /c = n — p if p G {0, 1, n — 1, n} or p{n — p) < m — (j{q) , ox ii C is, n generalized Gabidulin code 
or an ELS. 

Proof: The cases p G {0, n — l,n} are straightforward. In all other cases, since A; < n — p by 
Proposition [131 it suffices to prove that k > n — p. First, suppose p = 1, then k satisfies q'^^ > 
by the sphere covering bound. However, Vi{q"^,n) < q"^+"- < g^m^ hence k > n — 2. Second, if 
p(n — p) < m — a{q), then < ^ {p{n — p) + a{q)) < 1 and A; > n — p by Proposition [13] Third, if 
C is an {n,k,n — k + 1) generalized Gabidulin code with k < n, then there exists an {n, k + l,n — k) 
generalized Gabidulin code C such that C C C . We have p > df^{C') = n — k, as noted in Section Hl-CI 
and hence k > n — p. The case A; = n is straightforward. Finally, if C is an ELS of dimension k, then 
for all X with rank n and for any c G C, (iR(x, c) > rk(x) — rk(c) > n — k. ■ 
A similar argument can be used to bound the covering radius of the cartesian products of generalized 
Gabidulin codes. 
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Corollary 3: Let Q be an {n,k,d^) generalized Gabidulin code (n < m), and let be the code 
obtained by Z cartesian products of ^ for / > 1. Then the rank covering radius of C/' satisfies > ^r— 1. 

Note that when n = m, ^' is a maximal code, and hence Corollary |3] can be further strengthened. 

Corollary 4: Let Q be an (m, k, d^) generalized Gabidulin code over GF(g'"), and let be the code 
obtained by / cartesian products of Q. Then p{G^) = d^ — 1. 

The tightest bounds for the dimension of linear codes with given covering radius known so far are 
given in Table |ll] for (7 = 2, 4 < m < 8, 4 < n < m, and 2 < p < 6. 

E. Asymptotic covering properties 

Table |l] provides solutions to the sphere covering problem for only small values of m, n, and p. Next, 
we study the asymptotic covering properties when both block length and minimum rank distance go to 
infinity. As in Section [TVl we consider the case where lim„__»oo ^ = b, where 6 is a constant. In other 
words, these asymptotic covering properties provide insights on the covering properties of long rank 
metric codes over large fields. 

The asymptotic form of the bounds in ^ are given in the lemma below. 



Lemma 12: For < (5 < min{l, b'^}, v{5) = lim„_oo "^^"^''^ = 5(1 + 6 - b5). 

Proof: By LemmaQ we have < ^^^(g™, n) < qMm+n~dR)+a{q) _ taking the logarithm 

and dividing by n, this becomes 5{1 + 6 — 6(5) < loggm(yL5„j (g"^, n))/n < 6(1 + b — b5) + The 
proof is concluded by taking the limit when n tends to infinity. ■ 

Define r = ^ and k{r) = lim„_oo inf [ '"s.- The bounds in (HB and (ED together solve 



log ™ Vis^i {g"',n) 



the asymptotic sphere covering problem. 
Theorem 1: For all 6 and r, we have 

k{r) = {l-r){l-br). (23) 
Proof: By Lemma [T2l the sphere covering bound in ([TT]) asymptotically becomes k{r) > (1 — r)(l — 
br). Also, from the bound in (|2T]) . we have 

log K^iq"", n, p) < log ,„„ — — — - + 0((mn)"^ In(mn)). 

Fp(^™,n) 

By Lemma [T2I this asymptotically becomes k{r) < (1 — r)(l — 6r). Note that although we assume n < m 
above for convenience, both bounds in (fTTI) and ((2TI) hold for any values of m and n. ■ 



February 1, 2008 



DRAFT 



20 



VI. MacWilliams identity 

For all V G GF{q"^)"- with rank weight r, the rank weight function of v is defined as /r(v) = y^x"'~^. 
Let C be a code of length n over GF{q"^). Suppose there are Ai codewords in C with rank weight i 
(0 <i < n), then the rank weight enumerator of C, denoted as W(?(x,?/), is defined to be 

n 

W^{x, y) ^(^) = E ^iV'^'^''- 

vec i=o 

A. q-product of homogeneous polynomials 

Let a{x,y;m) = Yl\=Q(^i{^)y^^^~^ b{x,y\m) = Yl^j=Q^j{^)y^ be two homogeneous 
polynomials in x and y of degrees r and s respectively with coefficients ai{m) and hj{m) respectively. 
ai{m) and hj (m) for i, j > in turn are real functions of m, and are assumed to be zero unless otherwise 
specified. 

Definition 1 (q-product): The q-product of a{x, y; m) and h{x, y; m) is defined to be the homogeneous 
polynomial of degree (r + s) c{x, y; m) =^ a{x, y; m) * b{x, y; m) = X^^Iq Cu{m)y'^x'^~^'^~^, with 

u 

Cuim) = ^ q"-"" ai{m)hu-i{m - i). (24) 

1=0 

We shall denote the g-product by * henceforth. For n > the n-th g-power of a(x, y; m) is defined 
recursively: a{x,y;m)^^^ = 1 and a(x, y;m)["l = a(x, y; m)["~^l *a{x,y;m) for n>l. 

We provide some examples to illustrate the concept. It is easy to verify that x*y = yx, y*x = qyx, 
yx*x = qyx^, and yx * {q'^ — l)y = {q'"' — q)y^x. Note that x * y 7^ y * x. It is easy to verify that the 
g-product is neither commutative nor distributive in general. However, it is commutative and distributive 
in some special cases as described below. 

Lemma 13: Suppose a(x, y; m) = a is a constant independent from m, then a(x, y; m) * 6(x, y; m) = 
b{x,y;m) * a{x,y;m) = ab{x,y;m). Also, if deg[c(x, y; m)] = deg[a(x, y; m)], then [a(x,y;m) + 
c(x, y; m)] * b{x, y; m) = o(x, y; m) * 6(x, y; m) + c(x, y; m) * 6(x, y; m), and 6(x, y; m) * [a(x, y; m) + 
c(x, y; m)] = 6(x, y; m) * a(x, y; m) + 6(x, y; m) * c(x, y; m). 

The homogeneous polynomials a/(x, y;m) =^ [x + {q^ — l)y]''l and bi{x,y\m) =^ (x — y)''] are 
very important to our derivations below. The following lemma provides the analytical expressions of 
Oi(x, y; m) and bi{x, y; m). 
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Lemma 14: For i > 0, (jj == ^^^2^- / > 0, we have y^^ = q^'y'' and x''' = xK Furthermore, 

a{m,u)y''x^-'', (25) 



2 

/ 



ai{x,y;m) = ^ 

u=0 



bi{x,y;m) = ^ 



(26) 

Note that ai{x,y;m) is the rank weight enumerator of GF(g'")'. The proof of Lemma [141 is given in 
Appendix [Cl 

Definition 2 (q-transform): We define the q-transform of a{x,y;m) = Yll=o^ii''^)y^^^~^ '^he ho- 
mogeneous polynomial a(x, y; m) = X^^^g * a:''"^*^. 

Definition 3 (q-derivative [35]): For q > 2, the g-derivative at x / of a real-valued function /(x) 
is defined as 

def fjqx) - fix) 
{q-l)x 

For any real number a, [f(^x)+ag{x)]^^^ = f^^\x)+ag^^\x) for x / 0. For > 0, we shall denote the 
i/-th g-derivative (with respect to x) of f{x,y) as f^'^\x,y). The 0-th g-derivative of f{x,y) is defined 
to be /(x, 2/) itself. 

Lemma 15: For < z/ < (x')^'^) = (3{l^v)x^~^ . The i/-th g-derivative of f{x,y) = Yll=o fiy^^^~^ 
is given by /M(x,2/) = EI=o Also, 

a[''\x,y;m) = [3{l,v)ai^y{x,y;m) (27) 

bf\x,y-m) = (3{l,u)bi^^{x,y;m). (28) 
The proof of Lemma [15] is given in Appendix iDl 

Lemma 16 (Leibniz rule for the q-derivative): For two homogeneous polynomials /(x, y) = Y^\=o fiV^^^ 
and g{x,y) = Y^j^q gjy-' x'^~^ with degrees r and s respectively, the z^-th (u > 0) g-derivative of their 
(^-product is given by 



[f{x,y)*g{x,y)f^ = Y, 



1=0 

The proof of Lemma [16] is given in Appendix [E] Lemma [16] gives the z^-th g-derivative of g-products 
of homogeneous polynomials. 

The -derivative is similar to the g-derivative. 

Definition 4 (q~^ -derivative): For q > 2, the g-derivative at y / of a real-valued function g{y) is 
defined as 
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For > 0, we shall denote the z^-th q ^-derivative (with respect to y) of g{x,y) as g^'^^{x,y). The 
0-th -derivative of g{x, y) is defined to be g{x, y) itself. 

Lemma 17: For < < /, the z^-th g"^ -derivative of y^ is (y')^''^ = Also, 

al''\x,y;m) = l3{l,v)q~"" a{m,v)ai^^{x,y]m - v) (30) 

h\''\x,y-m) = (-ir/3(/,z.)6z_,(x,y;m). (31) 
The proof of Lemma [17] is given in Appendix IF] 

Lemma 18 (Leibniz rule for the q"^ -derivative): For two homogeneous polynomials /(x, y; m) = X][=o fiU^^^ 
and g{x,y;m) = Ylj=o djV'^ ^^'^ ^^^^ degrees r and s respectively, the z^-th (v > 0) (7^ ^-derivative of 
their g-product is given by 



[f{x, y- m) * g{x, y; m)]^"'' = ^ 



y. ^) ^ gi-n^a:, y-m- I). (32) 



1=0 

The proof of Lemma [18] can be found in Appendix [ 



B. The dual of a vector 

For all u, V € GY{q^)'^, let u - v denote the standard inner product of u and v. For any linear subspace 
C C GF(g™)", = {u G GF(g'")"|u • V = OVv G £} is referred to as the dual of C. 

As an important step toward our main result, we derive the rank weight enumerator of (v)^, where 



det 

V € GF(g™)" is an arbitrary vector and (v) = {av : a E GF(g'")}. Note that (v) can be viewed as 
an (n, 1) linear code over GY{q"^) with a generator matrix v. It is remarkable that the rank weight 
enumerator of (v)"*" depends on only the rank of v. 

Berger [36] has determined that the linear isometrics for the rank distance are given by the scalar 
multiplication by a non-zero element of G¥{q^), and multiplication on the right by an nonsingular 
matrix B G GF(g)"^". We say that two codes C and C are rank-equivalent if there exists a linear 
isometry / for the rank distance such that /(C) = C". 

Lemma 19: Suppose v has rank r > 1, Then C = (v)"*" is rank-equivalent to C x GF((/™)"~'', where 
C is an (r,r — 1,2) MRD code and x denotes cartesian product. 

Proof: We can express v as v = vB, where v = (wq, . . . , fr-i, . . . , 0) has rank r, and B G 
Qjr(g,)"X" has full rank. Remark that v is the parity-check of the code C x G¥{q^)'^~'^ , where C = 
((i^o, . . . , Vr-i))"*" is an (r, r — 1,2) MRD code. It can be easily checked that u € £ if and only if 
u = vlB^ G (v)-^. Therefore, (v)^ = CB^ , and hence C is rank-equivalent to (v)^ = C x GF(g™)"~^. 
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We hence derive the rank weight enumerator of an (r, r — 1, 2) MRD code. Note that the rank weight 
distribution of Unear Class-I MRD codes has been derived in [2], [8]. However, we will use our results 
to give an alternative derivation of the rank weight distribution of linear Class-I MRD codes later, and 
thus we shall not use the result in [2], [8] here. 

Proposition 15: Suppose G GF{q"^Y has rank r (0 < r < m). The rank weight enumerator of 
Cr = (v)"*" depends on only r and is given by 

Wl{x,y) = g-™ {[x + {q^ - 1)?/]^ + (g" - l){x - . (33) 

Proof: We first prove that the number of vectors with rank r in Cr, denoted as Ar^r, depends only 
on r and is given by 

Ar,r = q-"'Hm,r) + (q^ - (34) 

by induction on r (r > 1). Eq. (l34l ) clearly holds for r = 1. Suppose Eq. (l34b holds for r = f — 1. 

We consider all the vectors u = [uq, . . . ,Uf-i) G Cf such that the first f — 1 coordinates of u are 
linearly independent. Remark that Uf-i = —v^}^ Si=o completely determined by uq, . . . , Uf-2- 

Thus there are Nf-i{q"^, f — 1) = a{m, f — 1) such vectors u. Among these vectors, we will enumerate 
the vectors t whose last coordinate is a linear combination of the first f — 1 coordinates, i.e., t = 
{to, tr-2, Ya=o '^i^i) ^ whcrc Oj G GF(g) for < i < f - 2. 

Remark that t G >Cf if and only if (to, ■ • • , tf-2) ■ (^0 + CLQVf-i, ■ ■ ■ , Vf-2 + af-2Vf-i) = 0. It is easy 
to check that v(a) = (vq + ooff-i, • • • , + af_2t'f-i) has rank f — 1. Therefore, if ao, . . . , af_2 are 
fixed, then there are Af-i^f-i such vectors t. Also, suppose X][=o + ''^r-i YllZo biti = 0. Hence 
Yll=o(.^i ~ bi)ti = 0, which implies a = b since tj's are linearly independent. That is, (v(a))^ n 
(v(b))^ = {0} if a / b. We conclude that there are q^~^ Af-i^f-i vectors t. Therefore, Af^f = 
a(m,f - 1) - q^-^Af-i^f~i = g"'"[a(m,f) + (g*" - l){-lY q""^']. 

Denote the number of vectors with rank p in Cr as A^^p. We have Ar^p = [^j Ap^p [8], and hence Aj.^p = 
[a(m,p) + (g--l)(-l)Pg-^]. Thus, Wl^{x,y) = TJp=o^r,px''"'y'' = + (^^ " 1)^]^'^ + 

We comment that Proposition [TS] in fact provides the rank weight distribution of any (r, r — 1, 2) MRD 
code. 

Lemma 20: Let Co ^ GF(g'")^' be a linear code with rank weight enumerator Wq^ (x, y), and for s > 0, 

det 

let Wq {x,y) be the rank weight enumerator of Cs = Cq x GF{q''^Y. Then (x,y) only depends on 
s and is given by 

Wl{x,y) = Wl{x,y) * [x + {q^ - . (35) 
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Proof: For s > 0, denote W^S (x, y) = Xl^to We will prove that 



Bs,u = ^ Q^^Bo^i 



s 

u — i 



a{m — i,u — i) (36) 

7 A — /, 

i = 

by induction on s. Eq. (l36l ) clearly holds for s = 0. Now assume (l36l ) holds for s = s — 1. For any 
= {xq, . . . , Xr+s-i) S Cg, we define xj_i = (xq, • • • , Xr+j_2) € Then rk(xs) = n if and only if 
either rk(xs-i) = u and Xr+s-i G 6(xs-_i) or rk(xg_i) = n — 1 and x-r+j-i ^ S(xs_i). This implies 
Bs,u = q''Bs-i,u + (g'" - = q'~'Bo,i [,^] a(m -i,u-i). ■ 

Combining Lemma [T9j Proposition [151 and Lemma [20l the rank weight enumerator of (v) can be 
determined at last. 

Proposition 16: For v € GF(g'")" with rank r > 0, the rank weight enumerator of C = (v)"*" depends 
on only r, and is given by 

Wl{x,y) = q-"^{[x + {q^ - l)yp + {q^ - l)(x - y)M * [x + (g- - . (37) 

C. MacWilliams identity for the rank metric 

Using the results in Section IVI-BI we now derive the MacWilliams identity for rank metric codes. But 
first, we give two definitions from [25] that are needed in our derivation. 

Definition 5: Let C be the field of complex numbers. Let a € G¥{q^) and let {1, ai, . . . , a^.i} 
be a basis set of GF((7™). We thus have a = oq + aiai + . . . + am-iO-m-i, where ai E GF(g) for 

def 

< i < m — 1. Finally, let C £ C be a primitive g-th root of unity, x(o) = provides a mapping from 
GF(g™) to C. 

Definition 6 (Hadamard transform): For a mapping / from GF(g™)" to C, the Hadamard transform 
of /, denoted as /, is defined to be 

/(v)= E X(u-v)/(u), (38) 

uGGF(g"')" 

where u • v denotes the inner product of u and v. 

Let C be an {n,k) linear code over GF((7™), and let W(^(x,y) = Yll=o ^iV^ be its rank weight 
enumerator and VF^i(x,y) = Yl^=o BjU^'x"^^ be the rank weight enumerator of its dual code C"*-. 

Theorem 2: For any (n, k) Unear code C and its dual code over GF(g"^), 

W^c- y) = i^^^c + - 1)2/, ^ - 2/) > (39) 
where is the g-transform of W^. Equivalently, 

n n 

Bjy^x^-^ = q-^^ ^^(^ - y)^"^ * + - 1)2^]'""'' • (40) 

j=0 i=0 
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Proof: We have rk(Au) = rk(u) for all A G GF(g™)* and all u G GF(g™)". We want to determine 
/r(v) for all V G GF((7™)". By Definition [6l we can split the summation in Eq. (1381 ) into two parts: 

/r(v) = I;x(u-v)A(u)+ x(u-v)A(u), 
ue-C uGF\£ 

where C = (v)^. If u G £, then x(u-v) = 1 by Definition|5j and the first summation is equal to W^{x, y). 
For the second summation, we gather vectors into groups of the form {Aui}, where A G GF((7™)* and 
ui • V = 1. We remark that for u G F\C (see [25, Chapter 5, Lemma 9]) 

^ x(Aui.v)/,(Aui)=/,(ui) X(A) = -/r(ui). 

AGGF(g")* AeGF(g'")* 

Hence the second summation is equal to — qr}_i Wp^^{x, y). This leads to /r(v) = qr}_i [q^W^{x, y) — 
Wpix^y)]. Using Wp{x,y) = [x + {q^ — and Proposition [T6l we obtain /r(v) = {x — y)M * 

By [25, Chapter 5, Lemma 11], any mapping / from F to C satisfies J2v£C^ f(^) ~ ]^ SvgC /(^)- 
Applying this result to /r(v), we obtain (l39l ) and (l40l ). ■ 
Also, Bj's can be explicitly expressed in terms of Ai's. 
Corollary 5: We have 

1 " 

Bj = —Y,^iPjii-^m,n), (41) 

where 

(-l)'g'^'g'("-')a(m - /,i - /). (42) 

Proof: We have (x — * (x + (g*" — = ^"^q m, The result follows 

from Theorem |2l ■ 

Note that although (1411 ) is the same as that in [2, (3.14)], Pj{i;m,n) in (l42l ) are different from Pj{i) 
in [2, (AlO)] and their alternative forms in [37]. We can show that 

Proposition 17: Pj{x;m,n) in (l42l ) are the generalized Krawtchouk polynomials. 

The proof is given in Appendix |Hl Also, it was pointed out in [37] that p ^^.'^'^j is actually a basic 
hypergeometric function. 

Proposition [TT] shows that Pj{x; m, n) in (l42l ) are an alternative form for Pj{i) in [2, (AlO)], and hence 
our results in Corollary [5] are equivalent to those in [2, Theorem 3.3]. 



Pj{i;m,n) = ^ 



i 




n — i 


I 
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D. Moments of the rank distribution 

Next, we investigate the relationship between moments of the rank distribution of a linear code and 
those of its dual code. Our results parallel those in [25, p. 131]. 
Proposition 18: For < < n. 



E 

i=0 



n — I 

V 



A, 



3=0 



n- J 
n — V 



(43) 



Proof: First, apply Eq. ^ to C^. We obtain Ai = g^C^^") YTj=o BjPi{j; m, n), and hence 



1=0 



n — I 

V 



n — I 

V 



Pi{j;m,n) 



j=0 1=0 

We have Y.7=o [Vl ^^O'^ ^) = g™^""'') ["Z^] [38, (29)], and the result foUows. ■ 
Proof: First, applying Theorem |2] to C"*", we obtain 

n n 

j=o i=o 
Next, we apply the g-derivative with respect to x to Eq. (l44b times. By Lemma [15] the left hand side 

(LHS) becomes YJl=o Pi^i - i, u)Aiy^x'^-^-'' 



, while by Lemma [16] the right hand side (RHS) reduces 



to ^m(fc-n) Y]]=oBji^j{x,y), where 

V 

i^j{x,y) = [bj{x,y;m) *an-j{x,y;m)]^'''> =^ 



j('^-l)0-%f{x,y)*a^^Zl\x,y;m 



x,y;m) 



By Lemma[l5] h)\x,y;m) = (3{j,l)(x-y)'^ '1 and a[^''_^. ^(x, y; m) = I3(n- j,v -l)an~j-u+i 
It can be verified that for any homogeneous polynomial y; m) and for any s > 0, (6 * as)(l, 1; m) = 
q"^^b{l, l;m). Also, for x = y = 1, ^^^(1, l;™.) = We hence have 1) = for j > i/, 

and V, (1, 1) = [;] - J> - j)^™^"-'^) for j < u. Since /3(n - j> - j) = [^^zj] - J> - j) 

and I5{v,v) = - j> - j), V',(l, 1) = /3(z., Applying x = y = 1 to the 

LHS and rearranging both sides using fi{n — = u), we obtain ( l43] ). ■ 

Proposition [18] can be simplified if u is less than the minimum distance of the dual code. 

Corollary 6: Let be the minimum rank distance of C"*-. If < < d^, then 



E 



i=0 



n — I 

V 



m{k—u) 



Ai = q' 

Proof We have Bq = 1 and Bi = . . . = = 0. 
Using the g^^-derivative, we obtain another relationship. 
Proposition 19: For Q < v < n, 

V r .-1 

{-lyq^^aim - j,u - j)q^^''^^'^Bj. 



(45) 



i=u j=0 



n- J 
n — V 



(46) 
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The proof of Proposition [19] is similar to that of Proposition [TH and is given in Appendix Jl Similarly, 
when V is less than the minimum distance of the dual code, Proposition [19] can be simplified. 
Corollary 7: If < < d'^, then 

a{m,v). (47) 
Proof: We have Bq = I Bi = ■ ■ ■ = B^ = 0. ■ 
Following [25], we refer to the LHS of Eq. ( [43] ) and ( [46] ) as moments of the rank distribution of C. 

E. Relation to Delsarte's results 

Delsarte [2] also derived the Mac Williams identity for rank metric codes, and below we briefly relate 
our results to those by Delsarte. 

Delsarte [2] considered array codes with the rank metric. In [2], the inner product between two mxn 
matrices A and B over GF(g) is defined as A • B =* Tr(AB^). Two matrices A and B are orthogonal 
if x(A • B) = 1, where x is a nontrivial character of the additive group GF(g), and dual codes are 
defined using this orthogonality. Delsarte then established an analytical expression (cf. [2, Theorem 3.3]) 
between the rank distance enumerator of an array code and that of its dual. 

Clearly the definitions of dual codes are different in our work and [2]. However, we show below the 
two definitions collide when dual bases are used. With a slight abuse of notation, the inner products 
between two vectors and two matrices are both denoted by • and dual codes of both vector and array 
codes are denoted by ±. For all vectors x G GF((7™)", we expand x into a matrix with respect to 
the basis B = {Aj^o^ °^ GF(g™) over GF(g) and refer to the matrix {a;»j}™=o,o~^ ^B- That is, 
Xj = Yli^^ ^ijf^i for < j < n. For a code C of length n over GF(q^), we denote Cb == {x^ G 
GF(g)™^"|x € C}. Clearly, if C is an {n,k) linear code over GF(g'"), then Cb is an {mn,mk) linear 
array code over GF(g). 

Lemma 21: Let E = {ej^g^ and P = {^jj^To^ be dual bases of GF(g™) over GF(g). Then for any 
a, b € GF(g™)", Trace(a • b) = a^; • bp. 

Proof: We have a • b = ajbj = ^ ,■ CLijeib^jcj)!.. Applying the trace function on both 
sides, we obtain 

Trace(a -13) = ^^^ j6fcjTrace(ei0fe) = X] X] = Tr(a£;bp). 

j i k j i 

m 

Proposition 20: For an (n, k) code C over GF(g™) and dual bases E and P of GF{q"^) over GF((7), 



•m{k—v) 
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Proof: Let v G C^, then for any u G C, v • u = 0. Hence v^; • up = by Lemma |2T] and 
xi^E • up) = 1 for all up G Cp. Therefore v^; G {Cp)^ and {C^)e ^ {Cp)^. Since \{C^)e\ = \{Cp)^\, 
we conclude that {C^)e = {Cp)^. ■ 
Proposition [20] implies that our identity in Corollary |5] can be derived from Delsarte's identity in [2, 
Theorem 3.3]. Suppose C is an (n, k) linear code over GF(g"^) with rank weight distribution Ai and its 
dual code has rank distribution Let E and P are dual bases of GF(g™) over GF((7). Note that 
Cp and {C^)e have the same rank distribution as C and C"*-, respectively. Also by Proposition l20l C p and 
{C-^)e are dual array codes. Thus applying Delsarte's identity to Cp results in Corollary |5l Note that the 
rank distance enumerator and rank weight enumerator are the same for linear codes. 

Our results in this section are different from Delsarte's results in several aspects. First, Pj{i;m,n) 
in (|42l) are different from Pj{i) in [2, (AlO)] and their alternative forms in [37]. In Proposition [TtI we show 
that Pj{x;m,n) are actually the generalized Krawtchouk polynomials, and hence Pj{i;m,n) in (l42l ) are 
equivalent to Pj{i) in [2, (AlO)]. Second, our approach to proving the MacWilliams identity is different, 
and intermediate results of our proof offer interesting insights (see Lemma [T9l and Proposition [T6l). Third, 
our MacWilliams identity is also expressed in a polynomial form (Theorem O similar to that in [25], 
and the polynomial form allows us to derive other identities (see, for example. Propositions [18] and [T9] ) 
that relate the rank distribution of a linear code to that of its dual. 

VIL Conclusions 

In this paper, we investigate the packing, covering, and weight properties of rank metric codes. We 
show that MRD codes not only are optimal in the sense of the Singleton bound, but also provide the 
optimal solution to the sphere packing problem. We also derive bounds for the sphere covering problem 
and establish the asymptotic minimum code rate for a code with given relative covering radius. Finally, 
we establish identities that relate the rank weight distribution of a linear code to that of its dual code. 
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Appendix 

The proofs in this section use some well-known properties of Gaussian polynomials [32]: 



n 




n 




k 




n — 


K 


n 




n — 


1 


k 




k 





+ 



n — 1 
k 

' - 1 



+ 



qn-k _ I 

„n-k+l 



n 



n — 1 
k-l 
n — 1 
k-l 
1 



1 



n 
k-l 



n 




'k 




n 




n — I 


A 




I 




I 




n — k 



(48) 
(49) 
(50) 
(51) 
(52) 
(53) 

/v J L ' J L ' J L'"' ~ 

A. Proof of Proposition |S] 

We first establish a key lemma. 

Lemma 22: If z £ Z and < p < n, then 

\AnBi{z)\<Viiq"',n)-qP-' 
Proof: By definition of p, there exists c G C such that dji{z,c) < p. By Proposition |2l |i?i(z) n 
Bp-i{c) \ gets its minimal value for df_{z, c) = p, which is q^^^ [^] by Proposition |4l A vector at distance 
< p — 1 from any codeword does not belong to A. Therefore, Bi{z) H i?p_i(c) C Bi{z)\A, and hence 
\A^Bi{z)\ = \Bi{z)\ - \Bi{z)\A\ < Vi{q^,n) - |5i(z) n Bp_i(c)[. ■ 
For a code C with covering radius p and e > 1, 



(54) 



7 



def 



l)[\C\Vp{q^,n) 



-\C\Vp_^{q^\n)]-\ 

< e\A\ - {e - 1)\Z\ 

< e|A| - (e- l)|^nZ| 
= e\A\Z\ + \AnZ\, 

where §B follows from |Z| < |C|yp(g™, n) - g"^", given in Section lUlCl 



(55) 
(56) 



1 < Yl ^c(i?i(a))+ Yl Ec{Bi{a)) (57) 

aeA\Z aGAnZ 
aeA 
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where ([57]) follows from Lemma |9] and \An Z\ < EdAf] Z). 



aeAxeBi (a)nz 

= E E ^c.({x}) 

xeZae-Bi(x)nA 

= ^\AnB,{^)\Ec{{^}), 
xez 

where (1581 ) follows from the fact the second summation is over disjoint sets {x}. Using Lemma [22l we 
obtain 



7 < 



E (y^i^ 



™ n) - qP-^ 



^i(r,n)- 
Combining and ([551), we obtain (ITtI) . 



^c({x}) 

(lcl^p('?'",n) 



(59) 



B. Proof of Proposition [77] 

Proof: Given a radius p and a code C, denote by Pp{C) the set of vectors in GF(g'")" that are at 
distance > p from C. To simplify notations, Q == g"^" and Pp{C) == (5^^|Pp(C)|. Let us denote the set 
of all codes over GF(g'") of length n and cardinality K as 5x- Clearly \Sk\ = (^)- Let us calculate 
the average value of Pp{C) for all codes C G Sk'- 

h E ^.(^) = ^^-^ E 1^.(^)1 = ^^-^ E E 1 



\Sk\ 



' ^' xeFCe5/f|dR(x,C)>p 



1 



cgSk xeF|dR(x,c)>p 



L^-^E 

xeF 
Q-Vp{q"',n 
K 



Q-Vp{q^,n) 
K 

Q 

K 



(60) 
(61) 
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Eq. (l60l) comes from the fact that there are ^"j^ '"•*) codes with cardinality K that do not cover x. 
For all K, there exists a code C G Sk for which Pp{C') is no more than the average, that is: 



PpiC) < 



QY'(Q-Vpiq^,n: 
k) \ K 

It A ( „m 



Let us choose K 



< (l-g-V,(<7™,n)) 

i°gQ(i-Q-^v-,to"',n)) J + 1 SO *at KlogQ (1 - Q-^iq'^^n)) < -1 and hence 
Pp{C') = (1 - Q"^Vp(^™,n))^^ < g-i. It follows that |Pp(C")| < 1, and C has covering radius at 
most p. ■ 



C. Proof of Lemma U4\ 

The proof is by induction on I. Clearly all the claims hold for I = 0. Suppose y''"^' = q"^-^y''~^ for 
I > 1, then = yV~^^ *y = q''^^q°''-^y'' = q'^^yK The proof for x^^ is similar. 

Suppose Eq. (l25l ) holds for Z = T — 1. We have ai{x,y;m) = ai_i{x,y;in) * {x + {q^ — \)y) = 
EUaz-.y"^'""- By (El, we have 



7 - 1 

u 



a{m, u) + q 



T- 1 

u-l 



1) 



a(m, u — l){q' 



m—u+l 



1) 



7 - 1' 




"/ - 1" 




+ 


u-l 


u 





a{m, u) 



a(m, u), 



(62) 



where Eq. ([62]) follows Eq. (1491) . 

Suppose Eq. (l26l) holds for / = T— 1. We have y; m) = y; m)*{x — y) = EL=o huV^^^"^- 

By dm), we have 



T- 1 
T- 1 

u 



-U—lr 



-i)V"+9 



u-l 



+ 



l - 1 
U-l 



l - 1 
U - 1 

-lYq^- 



(63) 



where Eq. (1631 ) also follows Eq. (|49 
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D. Proof of Lemma [75] 

The proof is by induction on v. Clearly all the claims hold for u = {). The z^-th g-derivative of 
follows Definition |3] Suppose Eq. (|27] ) holds for v = v — 1, then 



af\x,y]m) = [13 - I) ai_y+i{x,y;m)] 



i=l 



l-V + l 
1 



l-v + l 

i 

l-v 

E 



j=0 



(1) 

a{m,l-v + l-i)y^-^+^-' 
a{m, l-v - i)y^~''~'-x' 



i 



(64) 



where Eq. ( [64l ) follows Eq. ( [53] ). 

Similarly, suppose Eq. (]28]) holds fox u = v — 1, hence 

(1) 



bf\x,y;m) = [/?(/, - 1) y; m)] 

= /3(/,P-l) 5] 



l-u 



i=0 



i=l 

7 - V 

I 



l-V+l 

i 



(^_-[^l-y+^-iq(Ti-ii+i-^yl-u+i-i 



J-i 



= P{l,v)bi-p{x,y;m). 

E. Proof of Lemma [76[ 

We consider homogeneous polynomials /(x, y; m) = X]I=o fiV^^^^^ ^(^) 2/5 = Z]i=o ^i?/*^;^"* 
of degree r as well as g{x,y;m) = J2j=o djV'' v{x,y;m) = J2j=o'^jy'' of degree s. First, 

we need a technical lemma. 

Lemma 23: If ti,. = 0, then 



1 u{x,y;m) 

— [u[x, y; m) * v[x, y; m)) = * v{x, y; m). 

X x 



If Vs = 0, then 



1 , , , , NX / N v(x, y; m) 

— (n(x, y; mj * y; m)) = uix, qy; m) * , 

X X 

Proof: Suppose Ur = 0, then = Y^^Zq Uiy^x'^'~^~^. Hence 



(65) 



(66) 



u{x,y;m) 



r+s-l / k 



* v[x, y;m) 



E E q''ui{m)vk.i{m - I) y^x'+'-^-^ 



k=0 \l=0 

— (u(x, y; m) * vix^ y; m)). 

X 
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Suppose Vs = 0, then = Yl^^'^jV'' ^ ^- Hence 

I \ r+s—l / k \ 

v[x,y;m) _ ^ ( J(s-i) J„. ;\ 1 ,,fc^r+s-i-fc 



u{x, qy; m) 



X 



fc=0 \«=o 



(m(x, y; m) * f (x, y; m)). 



We now give a proof of Lemma [161 

Proof: In order to simplify notations, we omit the dependence of the polynomials / and g on the 
parameter m. The proof is by induction on v. For = 1, we have 

1 



[f{x,y)*g{x,y)]^^'^ 



■ [f{qx, y) * g{qx, y) - f{qx, y) * g{x, y) 



{q - l)x 

+ f{qx,y) * g{x,y) - f{x,y) * g{x,y)] 



1 



{q - l)x 



[fiQX, y) * {g{qx, y) - g{x, y)) + {f{qx, y) - f{x, y)) * g{x, y)] 



f, N g{QX,y) - g{x,y) f{qx,y) - f{x,y) 
fiqx, qy) * -. + -. * g{x, y) 



{q - l)x 



{q - l)x 



" f{x,y) *g'^^\x,y) + *g{x,y), 



(67) 
(68) 



where Eq. (l67l ) follows Lemma |23] 

Now suppose Eq. (l29l ) is true for u = u.ln order to further simplify notations, we omit the dependence 
of the various polynomials in x and y. We have 



1=0 
v 



JP-l)(T-l) 



(1) 



1=0 
v 



^ q(i^-l){r-l) (^qr-l ^ g{P-l+l) ^ f{l+l) ^ g{u-l)\^ 



(69) 



u+l 



1=1 



+ q 



V 
I - 1 



1=0 

V 

= E 

i=\ 

= E 

where Eq. ^ follows Eq. ([68]), and Eq. dVOll follows Eq. (|49 



V 

l-\ 



P + 1 
/ 



(70) 
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F. Proof of Lemma |77| 

Proof: The proof is by induction on v. Clearly all the claims hold for v = {). The i/-th -derivative 
of follows from Definition ID Suppose Eq. (l30b holds for: u = v — \, then 



l-u+l 

/3G,i7-l)g--ia(m,i7-l) 

l-u+l 



1-1? + I 

i 



l-i) 
i - 1 



= P{l,i?)q''''-'a{m,u -1) ^ 
= P{l,r')q'^''a{m,r')ai-p{x,y;m - u). 
Suppose Eq. dSTT ) holds for v = v — 1, then 



i-p+i 



i=l 



l-V + l 

i - 1 



iyq'''x^-''+^~'q^-'y'-^ 



{-lYP{l,v)bi_p{x,y;m). 



G. Proof of Lemma [781 



We consider homogeneous polynomials /(x, y; m) = Yl\=o fiV^^^"^ ^^id ^(^^i 2/; = Z]i=o ^j?/*^;''"* 
of degree r as well as g{x,y;m) = J2j=o djV'' ^'^'^ v{x,y;m) = J2j=oVjy^ x^~^ of degree s. First, 
we need a technical lemma. 

Lemma 24: If uq = 0, then 



If vq = 0, then 



y; m) j * t'(x, y; ml j = (7 * v\x^ y;m — 1). 

y y 



1/ / ^ t w r \ v{x,y;m) 
— [u(x, y; m) * v(x, y; m)) = u(x, qy; m) * 

y y 



(71) 



(72) 
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Proof: Suppose uq = 0, then = ^[^^ Uj+ix'' ^ Hence 



A:=0 \l=Q 



r+s / k \ 

fc=l \Z=1 / 
y; m) * v{x, y; m)). 

y 



Suppose vq = 0, then '"^^'^'"^^ = Yl^j^o'^j+i^^ ^ ^U^- Hence 



y 

r+s-l / k 



" fc=0 \«=0 



r+s /k-1 \ 
fc=l \«=0 / 

-(u(x, y; m) * v(x, y; m)). 

y 



We now give a proof of Lemma [HI 

Proof: The proof is by induction on u. For u = 0, the result is trivial. For = 1, we have 

[fix, y, m) * g{x, y; m)]^^^ = [f{x, q'^y, m) * g{x, q'^y; m) - f{x, q~^y; m) * g{x, y; m) 

+ fix, q~^y; m) * g{x, y; m) - f{x, y; m) * g{x, y; m)] 
= ^^^i\-^^-^y [/(^' I'^y^ * (six, q~^y; m) - g{x, y; m)) 

+ {f{x, q~^y\ m) - f{x, y; m)) * g{x, y; m)] 
g{x,q-^y;m) - g{x,y;m) 



f{x,y;m) * 



(q-^ - i)y 



sf(.x,q ^y]m) - f{x,y;m) 

+ q 7 — ] 7T *g(x,y;m-l) (73) 

(g-i - l)y 

= f{x,y;rn)*g^^\x,y;m) +q''f^^\x,y;m)*g{x,y;m-l), (74) 
where (1731 ) comes from Lemma [24l 
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Now suppose Equation (I32b is true for u. In order to further simplify notations, we omit the dependence 
of the various polynomials in x and y. We have 

{1} 



[/(m)*5(m)]^^+i> = 



1=0 



.l{s-u+l) 



1=0 L -I 

V 

= E 

(=0 

+ E 



(75) 



V 



V 
I - 1 







V 






+ 


^ - 1 


) 





E 

1=0 



I 



(76) 



where (1751 ) comes from (1741) and (1761) comes from ([SU 



H. Proof of Proposition [77| 

Proof: It was shown in [37] that the generalized Krawtchouk polynomials are the only solutions to 
the recurrence 



+ l;m + + 1) = + 1; m, n) + 



m,n) 



(77) 
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with initial conditions Pj{0;m,n) = ["]a(m,j). Clearly, our polynomials satisfy these initial conditions. 
We hence show that Pj{i;m,n) satisfy the recurrence in Eq. (TTTl) . We have 

i+i r. , n r ^ . i 

(-l)'g'"'g'("-*)a(m + + 1-1) 



Pj+i(i + l;m + l,n + l) 



i + l 
I 



n — I 



'i + l 
I 



m + 1 — I 

J + 1-1 



E 

i+l 

E 

E{«' 

z=o 

(-l)'g"'g'("-*)a(n-i,i + l-/) 

i 

E 



Xj(f^^i^-i)a{n-i,2 + \-l) 



+ 



I - 1 



m — / 

.i + l - ^. 



+ 



m — / 



(78) 







m — / 






j + 1 - / 



1=0 

i+l 

+ E 

1=1 

i+l 

+ E 



i 


qj+i 


m — I 


I 







I - 1 



1=1 



I 

I - 1 



^J+l-l 

m — I 



(-l)'g'^'g'("-^)a(n-i,j + l-0 
-l)'g"'g'("-*)a(n-i,j + l-0 

iyq'^'q^i"-i)a{n-i,j + 1-1) 
(-l)V'g'("-)a(n-i,i + l-/), (79) 



m — / 
j + 1- I 



where (1781 ) follows from (I501 ). Let us denote the four summations in the right hand side of Eq. (1791) as 
A, B, C, and D respectively. We have A = q^'^^Pj^i{i;m,n), and 



« = E 
E 



C 



1=0 



m — / 



{-l)^q-^q^i^-i)a{n l){q''-'+^ - q^), 

iy+\'^'+^q(^+^)i^-i)a{n - i,j - I) 



(80) 



m — I — 1 



D 



1=0 
1=0 



m — I 
.j-l. 



i 




m — I 


I 




J-l. 



(-l)'g'"'+^g'(""*)Q(n-i,j-/)g 



^m—l 



1 



(81) 
(82) 
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where ([IB follows from dSB and ([82l) follows from both and dSB- Combining (HO]), (H), and (l82l) . 
we obtain 



B + C + D 



i=0 







m — I 


I 




.j-l. 




- - 



(-l)'g'^'+ig'("-*)a(n-i,j-/) 



tin— I 



-,m—l 



-q^ Pj{i; m, n). 



I. Proof of Proposition [79l 

Before proving Proposition [191 we need two technical lemmas. 
Lemma 25: For all m, v, and I, we have 



6{m,iJ,j) 



i=0 



{-iyq'''a{m-i,u) = a{u, j)a{m - j,u - j)q^^"'-^\ (83) 
Proof: The proof is by induction on j. When j = 0, the claim trivially holds. Let us suppose it 
holds for j. We have 

(-l)'(7°''a(m — i, v) 



5{m,v,j + 1) 



E 

1=0 

i+i 

E 

2 = 

3 

E 

i=0 

J 



i + 1 

i 



+ 



J 

i - 1 



{-lyq^^aim-i,!^) 
i+i 



(-l)^g"'a(m-i,i/) +^ 
j=i 

^q' ] {-iyq-^a{m-i,u)-j2 
1=0 ' 

j 

E« 



i=0 



J 

i - 1 



( — l)'(7'^'a(m — i, u) 



1=0 
^m— 1 /'„!^ 



-iyq'''a{m -\-i,v- \)q'^-^-\q^ - 1) 

(i- 'ist -\)b{m-\,v-\;2) 

= a{v~j + l)a(m - j-l,i/-J - 

where ([84] ) follows from (l50l ). 

Lemma 26: For all n, i/, and j, we have 



1=0 



'J 




'n - J 


I 




v-l 



n- J 
n — V 



(84) 



(85) 
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Proof: The proof goes by induction on j. When j = 0, the claim trivially holds. Let us suppose it 
holds for j. We have 

i+i 



1=0 

i+i 



1 + 1" 




"n - 1 - f 






u-l 



3 



/-I 



n - 1 - j 



i+i 







n - 1 - j" 






I/-/ 



«=1 



1 


J 




"n - 1 - j" 




/-I 







(87) 



where (186] ) follows from (|49l ). Let us denote the first and second summations in the right hand side of 
dSVl) as A and B, respectively. We have 

i 



'J 




"n - 1 - ] 


I 




u-l 



1=0 

n-l-j 



n — 1 — u 



(88) 



and 



d-i 



1=0 



n-l-j 
u-l-l 



1=0 



n-1- j 
v-l-l 



/(-"-)(-l)V'a(^-l-/,j-0 



= -q^+^-^ein-l.v-l-j) 
[ n — V 

Combining ([86l), ([88]l, and ([89l), we obtain 
where (l90ll follows from (l52l). 



(89) 



n-1 - / 
n — 1 — 
n - 1 - i 

n — 
n-l-j 
n — z/ 



J+n—u 



n-l-j 
n — V 



((7^-^-1)^-^ ^+9"-^^ 

gi^-j - 1 



(90) 
(91) 
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We now give a proof of Proposition \T9\ 

Proof: We apply the -derivative with respect to y to Eq. (l44l ) u times, and we apply x = y = 1. 
By Lemma [T7] the left hand side (LHS) becomes 

n I- .-1 



The right hand side (RHS) becomes q^^^-'^) Yl]=oBjil^j{l, 1), where 
i^jix,y) = [bj{x,y;m) * an-j{x,y]m)]^''^ 



(92) 



E 

i=0 



6j_i(x, y; m) * a(m - Z, i/ - l)an-j-u+i{x, y\m-v) 
hj-i{x, y\ m) * a(m -l,v - l)an~j-u+i{x, y;m- v), 



(93) 
(94) 



where (1931 ) and (l94l) follow from Lemmas [18] and [T7] respectively. 
We have 



i 



E 

u=0 



j - I' 

i 



{-lyq^'-aim-i-l,!^ -I) 
{-iyq"'a{m-l-i,i^-l) 



^ qi{n-j-u+l) 
i=0 

(m—v){n—v—j+l) ^ ^ 



n — j — V + I 
u — i 



a{m — v — i,u — i) 



(95) 



where (1951) follows from Lemma |25] Hence 



E 







"n - i" 









n- J 
n — V 



(96) 



where ( |96l ) follows from Lemma l26l Incorporating this expression for ipj {1,1) in the definition of the 
RHS and rearranging both sides, we obtain the result. ■ 
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m 


n 


p=l 


P = 2 


p = 3 


p = 4 


p = 5 


p — 6 


2 


2 


b 3-4 A 


1 










3 


2 


b 4 B 


1 












3 


b 11-32 C 


a 2-4 C 


1 








4 


2 


b 7-8 B 


1 












3 


b 40-64 B 


b 4-8 C 


1 










4 


c 293-1024 C 


b 10-64 C 


a 2-8 C 


1 






5 


2 


b 12-16 B 


1 












3 


b 154-256 B 


b 6-8 B 


1 










4 


b 2267-4096 B 


b 33-256 C 


a 3-8 C 


1 








5 


b 34894-21'^ C 


b 233-2979 E 


b 10-128 C 


a 2-8 C 


1 




6 


2 


b 23-32 B 


1 












3 


b 601-1024 B 


a 10-16 B 


1 










4 


b 17822-2^^ B 


b 123-256 B 


b 6-16 C 


1 








5 


b 550395-2^" B 


b 1770-2" C 


c 31-256 C 


a 3-16 C 


1 






6 


c 17318410-2^6 C 


c 27065-424990 E 


c 214-4299 E 


c 9-181 D 


a 2-16 C 


1 


7 


2 


b 44-64 B 


1 












3 


b 2372-4096 B 


a 19-32 B 


1 










4 


b 141231-2^** B 


c 484-1024 B 


b 10-16 B 


1 








5 


b 8735289-2^^ B 


b 13835-21'^ B 


b 112-1024 C 


a 5-16 C 


1 






6 


b 549829402-2^° B 


c 42229-2^^ C 


b 1584-2^^ C 


b 31-746 E 


a 3-16 C 


1 




7 


b 34901004402-2^^ C 


c 13205450-244855533 E 


b 23978-596534 E 


c 203-5890 E 


a 8-242 D 


a 2-16 C 



TABLE I 



Bounds on KRlq"", n, p), for 2<m<7, 2<n<m, AND 1 < /9 < 6. FOR EACH SET OF PARAMETERS, THE TIGHTEST 
LOWER AND UPPER BOUNDS ON KR{q"^,n, p) ARE GIVEN, AND LETTERS ASSOCIATED WITH THE NUMBERS ARE USED TO 
INDICATE THE TIGHTEST BOUND. THE LOWER CASE LETTERS a-C CORRESPOND TO THE LOWER BOUNDS IN (fTTT ). (fTSl) . 
AND ([TT] ) respectively. THE UPPER CASE LETTERS A-E DENOTE THE UPPER BOUNDS IN ifTTT ). ( fTSl ). (fT9] |. (l20t . AND (1211 

RESPECTIVELY. 
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m 


n 


p = 2 


p = 3 


p = 4 


p = 5 


p = 6 


4 


4 


1-2 


1 









5 


4 


1-2 


1 











5 


2-3 


1-2 


1 







6 


4 


2 


1 











5 


2-3 


1-2 


1 









6 


3-4 


2-3 


1-2 


1 





7 


4 


2 


1 











5 


2-3 


1-2 


1 









6 


3-4 


2-3 


1-2 


1 







7 


4-5 


3-4 


2-3 


1-2 


1 


8 


4 


2 


1 











5 


3 


2 


1 









6 


3-4 


2-3 


1-2 


1 







7 


4-5 


3-4 


2-3 


1-2 


1 




8 


5-6 


3-5 


2-4 


1-3 


1-2 



TABLE n 

Bounds on k for g = 2, 4<m<8, 4<n<m, and 2 < p < 6. 
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